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Abstract 

The planar motion of a testing body in the field of gravity of the massive ring is studied. The 
"perihelion" precession due to non-relativistic reasons is demonstrated on this simple example 
of planar motion in the radially symmetric (in the plane of motion) situation. The ring as the 
source of the gravitational field can serve as the toy model for quasi-planar sources of gravity 
appearing in the galactic dynamics, etc. The results are used to estimate the real perihelion 
shift of the planet Mercury due to gravitational quadrupole moment of the Sun. 

1 The potential 

Let us consider the thin (one dimensional) homogenous massive ring of radius R with total mass 
M . We consider the cartesian coordinates x, y, z such that the plane Oxy coincides with the 
plane defined by the ring itself, and the origin of our coordinates coincides with the center of the 
ring. Our first task is to compute the Newtonian potential V = V{x,y) of the gravitational field 
of the ring in the plane Oxy. Because of the rotational invariance of the problem, the potential 
is the function of radial distance r = y^c^^T^ only. According to the Newton's gravitational 
law we have 
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We see that the potential has the singularity at r = i? (at the ring). It holds: 
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and tne sign oi equality takes place only tor r = n. We remind tne dennition ot tne complete 
elliptic integral of the first kind 
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with help of which we can express our potential as follows 
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Finally, we obtain the potential of the ring in the form 
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Having the potential we can derive the gravitational force that acts upon a testing body of mass 
m: 



f{r) = mg{r) = —mV'{r) 
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where E is the complete elliptic integral defined for < x < 1 by 

E{x) = I \l 1 — X siv? {y)dy 
Jo 

and we have denoted by g the gravitational acceleration. 




Figure 1: Plot of the potential (jT} as a function of the fraction r/R with 2nM / {nR) = 1. 

For the next purposes we will write down the approximative expressions of the potential in 
the following three cases: 
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where 75; « 0.577216 is the Euler's 7 constant and F is the Euler's Gamma function. 

Equation © tells us that at large distances from the ring the field is almost identical with the 
field of point mass M and the leading correction term (proportional to {r/R)"^) is the projection 
of the quadrupole field into the plane of the ring. Equation (@J shows that the center of the ring 
is an unstable equilibrium. If a testing particle is displaced a little bit from the center then it 
is forced to move towards the ring, and the force is proportional to the distance from the center 
(inverted linear harmonic potential). Finally, equation ^ shows the feature of the potential 
close the ring. Namely, there is the logarithmic divergency of the potential as r ^ i?. 



2 Equations of motion and the conservation laws 

An easy way how to obtain equations of motion is to write down the Lagrange function (the 
Lagrangian) of the testing particle with the mass m in polar coordinates (r, (/>) (we will consider 
the motion in the plane of the ring only, as we have already mentioned above): 

L = —m (r^ + r^(^^'] — mV{r), (6) 



where the dot stands for the time derivative, and derive the equations of motion as the Euler - 
Lagrange equations 

r + V'{r) = 0, r'^4> + 2rr<j) = 0. (7) 

The order of this system of two differential equations can be reduced using the conservation 
laws. First of all, the energy E is conserved: 

(^r^ + r^c/)^^ + mV{r) = E. (8) 

The angular momentum vector must have the form (0,0,i2). Since L does not depend on 0, we 
have the following expression for the z- component of the angular momentum: 

= ^- (9) 

The previous two equations can be combined into one expressing the radial velocity r as the 
function of the radial position r and the conserved quantities E and C: 
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ihe inequality 

— -2nr)-4^>0 

defines, for given values of E and £, classically accessible region, i.e. the minimal and the 
maximal tm value of radial coordinate of considered particle. 



3 Radial motion (free fall) 

The radial motion or the free fall is characterized by the condition £ = 0. We will consider the 
situation when the initial (radial) velocity equals zero. This means that E/m = F(ro), where 
ro is the initial position, and equation (|Tn|l becomes 

\r\=^[y^rQ)-V{r)f\ 

Our task is to compute the time of the fall to the ring T. First, we will consider the case when 
ro > R. Thus, 

V2Jr [Viro)-V{r)\ 

Explicitly, 

T(r„) = 1 



If the initial position of the particle is inside the ring (ro < R) then 

thus (|THl holds also in this case. The free fall at the ring from outside can be compared with the 
radial fall at the point mass M placed in the center of the ring. The time of duration of such a 
fall is given by 
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The quantities T and Tp are compared in the fig. |21 
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4 Circular orbits 

For a testing body on a circular orbit the centripetal force fc must be equal to the gravitational 
force Q. This means that 

fc = mac = -m(t?r = --^ — ^ = g. 

The last equation can obviously be fulfilled for any r > R but, of course, cannot be fulfilled for 
any r < R, i.e. the testing particle can orbit the ring along any circle from outside. 
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Figure 2: Left graph shows the dependence of the free fall duration T on initial distance tq/R 
according (fTT| . where wee have chosen for definiteness: kM = 1 and R = 1. Right graph shows 
the fraction T/Tp with the parameters chosen in the same way as in the left graph. We see that 
the quantity T approaches relatively Tp as tq/R grows. More detailed analysis would show that the 
difference Tp — T approaches a positive constant of the order of unity as tq/R runs to infinity. 



5 Small deformations of circular orbits - "perihelion" 
shift 

In this section we will consider the close-to-circular trajectory of a testing particle sufficiently 
distanced from the ring. This means that the potential of the ring can be approximated by the 
first two terms of (jSI): 

kM ( 1 kM 7 1 9 

y~ 1 + = + 4, 1 = --kMR^. 13 

It is well-known that the closed non-circular trajectories of testing particle in a central field exist 
only in two cases: F ~ and V ~ — 1/?", see eg. [J or detailed proof in (2]. Therefore our 
close-to-circular orbit will not be closed, however for a small deviation from the circle it can be 
regarded as slowly rotating ellipse. We are interested in the angle connecting the two nearest 
positions of the particle and the ring (the point of closest approach of a planet (from our solar 
system, of course) and the Sun is called perihelion) and the center of our coordinate system. 
Combining the eqs. © and (|Tn|l we easily can derive that the change of the polar angle Ac;/) 
between the two forthcoming closest approaches of the particle with respect to the ring is given 
by 



-2^/ ^^rn{E-V)-^dr, (14) 

where r™, and tm are the minimal and the maximal values of the radius vector of the particle, 
respectively. The angle difference IS.(f) can be written as: A(/) = 27r + where the difference 
27r comes from the Newton potential and additional term dcj) is called perihelion shift in the 
mechanics of planetary motion around the Sun or the periastron shift in the mechanics of binary 
stars. An illustration of the "perihelion shift" in our case is given in fig. |31 



5 



Y 



y 




Figure 3: Plot of the trajectory of a testing particle in the field of massive ring. The radius of the 
ring is -R = 1 and E and C were chosen so that the related Keplerian trajectory would be described 
by = 3 and vm = 15. The body of the ring is depicted as the central circle. The particle orbits the 
ring anti-clockwise, we see, from the left figure, that the "perihelion" moves in the same direction. 



Using (|T^ and the Taylor expansion of the subintegral expression in ifTHl we obtain 

where r = r(0) is the elliptic trajectory in the Newtonian field of a point mass. For this trajectory 
we have 
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where e is the eccentricity and p is the parameter of the ellipse with the semi-axis 
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Finally, we have 
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or using the relation 
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we can express the perihelion shift H16|l with help of geometrical parameters of considered per- 
turbed ellipse: 
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We have studied in details the motion of a testing particle in the gravitational field of the massive 
planar ring. As the main aim, we have determined the shift of the perihelion of the testing particle 
with non-circular orbit with sufficiently large radius. For r/R^l the gravitational field of the 
ring is practically identical with the field of central monopole of the mass M and the quadrupole 
field, see (0). The formula lfTH|l describes, in fact, the perihelion shift for the particle in the 
planar motion due to quadrupole perturbation of Newton potential. Hence, it can be uses to 
estimate the perihelion shift of the planet Mercury due to Solar quadrupole moment, considering 
the approximation that the trajectory of Mercury is placed in the plane of the solar equator. 
Following the results about the multipole expansion of the gravitational field of the Sun we have 



^0 = |l - J2^P2{com) + O {Rl/r 



where Mq is the mass of the Sun, Rq is the mean radius of the Sun, P2{z) = — l/2-|-3/2z^ is the 
Legendre polynomial that depends on the azimutal angle and J2 = (2.18 ±0.06) x 10^'' is the 
dimensionless quadrupole moment of the Sun - for a long time the value of solar gravitational 
quadrupole was very uncertain, this value is from jH]. In the equator of the Sun plane we have 
= tt/2 and therefore in the mentioned plane the potential has the form 



So, we can use the formula ifTHjl to compute the Mercury's perihelion shift due to the Sun's 
quadrupole moment identifying the constant 7 as follows: 

7 = -^J^kMqRI. (17) 

Finally, taking into account the Mercury's parameters: the mass m ^ 3.302 X W^^kg, distance 
from the Sun in the perihelion r^, ~ 4.6 x 10''/cm, orbital velocity in the perihelion vm ~ 59km/ s 
{C = mrmVM, orbital period is about 88 days), and the mass of the Sun: Mq « 2 x l{)^^kg and 
its radius Rq 7 x lO^fcm, we have 



3vrJ2 Q Q = 3^J2 ■ ® ® pa 3.3 X 10^^°. 18 



This means that during the period of 100 years, the perihelion of the Mercury is shifted in the 
direction of planetary rotation around the Sun in approximately 0.28 arcsec. This value can be 
compared with the well-known value of the relativistic shift of the perihelion that is about 35 
arcsec per a century. Our value (|T8|l is quite less that is in accordance with accepted fact that 
the Solar quadrupole does not effect essentially the Mercury's orbit in the time interval about 
100 years. 
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